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ABSTRACT 

Combined with other CMB experiments, the WMAP survey provides an ac- 
curate estimate of the baryon density of the Universe. In the framework of the 
standard Big Bang Nucleosynthesis (BBN), such a baryon density leads to pre- 
dictions for the primordial abundances of ^He and D in good agreement with 
observations. However, it also leads to a significant discrepancy between the 
predicted and observed primordial abundance of ''Li. Such a discrepancy is of- 
ten termed as 'the lithium problem'. In this paper, we analyze this problem 
in the framework of scalar-tensor theories of gravity. It is shown that an ex- 
pansion of the Universe slightly slower than in General Relativity before BBN, 
but faster during BBN, solves the lithium problem and leads to ^He and D pri- 
mordial abundances consistent with the observational constraints. This kind of 
behavior is obtained in numerous scalar-tensor models, both with and without 
a self-interaction potential for the scalar field. In models with a self-interacting 
scalar field, the convergence towards General Relativity is ensured without any 
condition, thanks to an attraction mechanism which starts to work during the 
radiation-dominated epoch. 
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Introduction 



The Big Bang Nucleosynthesis prediction for the "^He primordial abundance is tradi- 
tionally considered as a good estimate of the baryon density of the Universe. However, the 
recent measurement of the Cosmic Microwave Background (CMB) anisotropies by WMAP 
( iBennett fc all 120031 ) now provides, when combined with other CMB experiments (CBI and 
ACBAR), another accurate estimate of the baryon density: Qh ff = 0.0224 ± 0.0 009 (or 
1] X 10^° = 6.14 ± 0.25 in terms of the baryon to photon ratio) ( Spergel fc al 20031 ). This 



value of rj can be used to compute the primordial abundances of light elements (mainly ^He, 
■^He, D and ^Li). The comparison of predicted and observed abundances is then a way to 
test the concordance between BBN and the CMB data. 

When one assumes the WMAP estimate of the baryon to photon ratio, the primordial 
abundances predicted by the standard BBN ( i.e. three families of light neutrinos, a neutron 
mean life r„ = 885.7 ± 0.8s, gravitation describ ed by General Relativity and a homogeneous 
and isotropic Universe) are (ICyburt &: alll2003l ): 



D/H 
^Li/H 



0.2484: 



h0.0004 
-0.0005 



2.75- 



-0.24 
-0.19 



X 10" 



while the observed abundances are: 



Yp 
D/H 
^Li/H 



0.2391 ±0.0020 
0.2452 ±0.0015 

2.78trst X 10-^ 
1.23l^-^^ X 10 
2.19 



ref. fLuridiana fc al 20031 
ref. fizotov fc al 19991 



Q3g xw ref. (iKirkman fc alll2003l ) 



0.68 ^ 1 n-10 
0.32 

0.46 in-lO 
0.38 ^ 



ref. fRvan fc al 2000) 
ref. f Bonifacio fc al 2002) . 



Here, the two va lues fo r the observed ¥„ correspond to independent estimates in refs. 
( iLuridiana &: all l2003l ) and (llzotov &: all Il999l ) based on observations of metal-poor extra- 
galactic ionized hydrogen regions. The D abundance is determined by observations of re- 
mote cosm ological clouds on the line of sight of high redshift quasars. Finally the ^Li/H 
estimate of iRvan fc all (120001 ) has been performed by observing halo stars, while the value of 
Bonifacio fc al ( 20021 ) comes from the observation of sta rs in the globular cluster NGC 6397. 
Although some recent estimates (IBoesgaard fc alll2005l ) leads to a somewhat smaller value. 
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sever al other independent determinations (iRyan fc a]||l999l : lThevenin fc a]||200ll : lAsplung &: al 
20051) obta i n obs erv ed primordial a b unda nces of ^Li/H similar to those reported either by 
Rvan fc all J2000h or ISonifacio fc all J2002h . 



We can see a relatively good agreement for the predicted and observed abundances of ^He 
and D, but a large discrepancy for ^Li/H. There had been many attempts to account for the 
low ^Li abundance indicated by observations. The first and most conservative possibility is 
the existence of systematic uncertainties in the observational determination of th e ''Li abun- 
dance . Ho wever, such unce r tainti es are not large enough to solve the problem (see lRyan fc al 
fl2000h and ISonifacio fc all toO'^ ). [n the same way , the s ystematic study of uncertainties 
on the nuclear reaction rates performed in ICoc fc all (120041 ) indicates that these uncertain- 
ties cannot explain the large discrepancy emphasized before. Modifications of the standard 
nucleosynthesis scenario then arise as possible ways to reconcile the predictions with the 
observations. Inhomogeneous nucleosynthesis has been analyzed (jJedamzik fc Rehnj|200ll ) 
but overproduces ''Li. Late particles decays can deplete ^Li, but they cannot account for the 
obser ved D and D/^He primordial abundances (jJedamziklll2004j : |Jedamzik2ll2004l : lEUis fc al 
20051 ) : another possibility is provided by the creation of baryons after BBN, accompanied 
by a lepton asymmetry before BBN, the two processes arising from Q-balls (jlchikawa fc al 
2004h . 



In this work, we will a ddress the lithium abundance problem i n the framewor k of scalar- 



tenso r theories of gravity flBrans &: Dickd Il96ll : iBergmannl Il968l : IWagonerl Il970l : iNordvedt 



19701 ). In these theories, gravitation is modified by the introduction of a scalar degree of 
freedom, which does not affect the standard nucle ar and particle phy sics. Such modifications 
arise as low-energy limits of superstrings theories (IGreen fc allll988l ) and they provide a way 
to change the expansion history of the Universe with minimal assumptions. We will show 
that such theories contain a mechanism that could be responsible for a low lithium abun- 
dance, despite the high baryon to photon ratio implied by the WMAP estimate. Big Bang 



Nucleosy nthesis in the context of a sca - 
the past (Wagoner 1973 : Barrowl 



1978 



ar-tensor cosraology has been extensively studied in 



Serna &: Alimi 1996: Alimi &: Sena 1997: Serna &: al 



20021 : IClifton fc all l2005l : Coc &: al 2006 ) , but the main goal of these works was to constrain 
the parameters of the theory thanks to the primordial abundances of ^He, D and ^Li, or 
to obtain the observed abundances with a matter density of the Universe different from its 
commonly assumed value. Moreover, general self-interacting scalar fields were not considered 
and, as we will see later, the introduction of such terms in the lagrangian can provide the 
mechanism necessary to explain the low ^Li abundance. 

The paper is arranged as follows. In Section [2l the scalar-tensor gravity theories and 
the implied cosmological models are presented. Then, in Section [3l we analyze the lithium 
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problem in the framework of different kinds of scalar tensor theories. We show that the 
solution of this problem requires a non trivial dynamics for the scalar field at the epoch of 
BBN. Finally, we discuss the generality of that solution. 



2. Scalar-tensor cosmological model 

2.1. Equations and observable quantities 

In scalar-tensor theories of gravity, the dynamics of the Universe contains a new scalar 
degree of freedom that couples explicitly to the energy content of the Universe. In units of 
c = 1, the action generically writes, in the Einstein frame: 

1 f f R* 1 



+ Srr,{^m,A\ip)gp, (2) 

being a bare gravitational constant, (p the scalar field, V{(f) its self-interaction term and 
A{ip) its couphng to matter. The functional Sm{'^m, ■A^{'^)g%) stands for the action of any 
field that contributes to the energy content of the Universe. It expresses the fact that 
all these fields couple universally to a conformal metric g^^ = A'^{(fi)g^j^, then implying that 
the weak equivalence principle (local universality of free fall for non-gravitationally bound 
objects) holds in this class of theories. The metric gfj,^ defines the Dicke- Jordan frame, in 
which standard rods and clocks can be used to make measurements (since in this frame, the 
matter part of the action acquires its standard form) . Defining 

and considering the transformations: 



|3 + 2cu(0)| = a-\(p) 



the action in the Dicke- Jordan frame writes: 



+ Smii^rn^Qtiu) (4) 
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Despite the conformal relation, these two frames have a different status: in the Dicke- 
Jordan frame, where the gravitational degrees of freedom are mixed, the lagrangian for the 
matter fields does not contain explicitly the new scalar field: the non gravitational physics has 
then its standard form. In the Einstein frame, the scalar degree of freedom explicitly couples 
to the matter fields, then leading for example to the variation of the inertial masses of point- 
like particles. Of course, the two frames describe the same physical world. Nevertheless, 
the usual interpretation of the observable quantities is profoundly modified in the Einstein 
frame, whereas it holds in the Dicke- Jordan frame, where the rods and clocks made with 
matter are not affected by the presence of the scalar field. That is why we will refer to the 
Dicke- Jordan frame as the observable one. 

Varying the Einstein frame action with respect to the fields yields the equations: 

Rl.-lR^'g;. = 87rG.T;, + T;, (5) 

□ V = -4'kGMv)T. + ^^^^ (6) 
V:t:^ = a{^)nv;^ (7) 

where T^, is the trace of the energy- momentum tensor of matter fields T^,^j,, related to the 
observable one by T^,^ = A-'^{(f)T^f,^, and T^f^ = 2(f^^Lp^^ - g*^^ig'^^'fa'fi3) - 2V{(f)g*^^ is the 
energy-momentum tensor of the scalar field. It is important to note that these equations 
reduces to those of General Relativity in presence of a scalar field if a{(p) = 0. 

Except when the contrary is explicitly stated, all the computations presented in this 
paper have been performed in the Einstein frame, because it leads to well posed Cauchy 
problems (that is elliptic and/or hyperbolic equations with a set of initial conditions) and 
has a perfectly regular dynamics. Nevertheless, the cosmological evolution resulting from 
these computations was later expressed in the Dicke- Jordan frame, where the interpretation 
of the observable quantities is easier. 



2.2. Homogeneous and isotropic Universe 

Under the assumption of a Universe filled with various homogeneous and isotropic matter 
fluids, the metric reduces to the Friedman-Lemaitre-Robertson- Walker (FLRW) form in the 
observable Dicke- Jordan frame: 

1 — kr^ 
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and, also in the Einstein frame because of the conformal relation: 

dsl = -dtl + al{t*)df 

with dt = A{ip(t^:))dt:^ and a(t) = A{ip{t^))a^:{t^:). 

Then, defining = (l/a*)((ia^,/(it*), the fields equations ([5]) become: 



2 8nG^:p^: 1.2 2 k 



-3— = AttG,{p, + 3p,) + 2ip^ - 2V{ip) (9) 



a* 



^ + 3H,ip+'^^^ = -ATTGM^)ip*-^P*) (10) 
dip 

where dots denote derivatives with respect to t*, whereas and are the total mass-energy 
density and pressure, respectively. They are conformally related to the energy density p and 
the pressure p in the Jordan frame by = A'^{ip)p and = A'^{(f)p. The observable Hubble 
parameter H{t) is related to the Einstein frame quantities by 

H{t) = -^(HM + ai^iQMQ) (11) 



2.3. Dynamics of the scalar field 

In the form given above, the time evolution of the scalar field is coupled, both in the 



Dicke- Jordan and the Einstein fr ame, to that of the scale factor. Previous works (jSerna fc al 



20021 : iDamour fc Nordvedtlll993l ) have found that, by introducing an appropriate change of 
variables, it is possible to obtain an evolution equation for the scalar field which is inde- 
pendent of the cosmic scale factor. Such an equation allows for the dynamical analysis of 
any scalar tensor theory and, in particular, its asymptotic behavior at early and late times. 
When the latter behavior implies a convergence mechanism towards General Relativity, one 
can ensures that the Solar System constraints will be satisfied. We will now extend such 
previous works to the general case of scalar-tensor theories with a self-interaction term. 

Introducing a new evolution parameter A = Ina^. + const, so that dX = H^dt^,, and 
denoting u' = du/dX, the evolution equations fl5])- ffTU]) lead to: 

2(1 — e + r?) ,, , ,, , , 

S — (f'^ 

= -e - (1 - 3wb)a{^) (12) 
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where 

e(A) 
r/((^, A) 

Wb 

e((^,A) 



3A; 



8nG,p,iX)al{X) 
p* 

dip 



4ttG^P^{X) 



This equation is similar, but not equivalent, to the motion equation of a mechanical 
oscillator with a varying "mass" rrieff = 2(1 — e + ri)/{3 — ip'"^), a varying "friction" i/g// = 
{1 — Wb — 4e/3 + 21]) and a "force" term F^ff = — O — (1 — 3wb)a. Then, one can define an 
effective potential Veff{(p,X) that verifies the relation: 

^^^^^^ = e(v.,A) + (l-3^,)a(^) (13) 

In the following, we will only consider fiat cosmologies (e = 0). In these cases, the 
"friction" term {1 — Wb + 2,7]) is always positive, and the dynamics of the scalar field is 
then analogous to a damped oscillating motion in the effective potential Ve//(v?, A). Such an 
effective potential (see Eq. f[T^ ) presents two different parts: a term due to the coupling 
function a{ip) and another term due to the self-interaction V{ip) of the scalar field. 

During the matter-dominated era {wb = 0), both terms are important to determine the 
effective potential. The minima of Vg// will be attractors for the dynamics of the scalar 
field and determine the behavior of the theory at late times. Consequently, if the relativistic 
value <y9 = is a minimum of this effective potential, the theory will converge towards General 
Relativity. 

Instead, during the radiation-dominated era [wb = 1/3), the only non- vanishing contri- 
bution to the effective potential is that due the self-interaction of the scalar field. Equation 
ffT^ then reduces to: 

|^¥^" + ^(l + 3r,)^' = -e (14) 

At very early times, when p^ is very high, B ~ 77 ~ so that self-interaction has a 
negligible effect on the scalar field dynamics. In this case (0 ~ ~ 0), the integration of 
equation IHM gives: 

2 
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where /c is a constant related to the value of ty^' at A = through /c^ = v'q /(3 — (Pq). Note 
that the above equation implies that — < ip' < \/?>. 

One can see from equation f[T^ that the velocity of the scalar field exponentially de- 
creases with A at very early times. When the initial (A = 0) velocity is not very high (i.e., 
when it is not close to -\/3), one can assume that the scalar field reaches the epochs of interest 
(those prior to BBN processes) with an almost vanishing velocity: 

= (before BBN). (16) 

For the sake of simplicity, we will consider throughout this paper a vanishing initial value 
of the scalar field velocity (see the Appendix for some examples where such a condition has 
been relaxed). Hereafter, all quantities at this initial time (prior to BBN) will be expressed 
with a subscript 'init'. 

The condition (flGl) does not however imply a constant (p value during the whole radiation- 
dominated epoch. Since p^. is a decreasing function of A, the contribution of V{ip>) to the 
effective potential will become non- negligible at some 'time' A. If V{(p) is chosen to be a 
function having a minimum at <^m; the scalar field will start to be attracted towards ipm- 
For example, the simple choice of a power law V{ip>) oc ^9^" for the self-interaction term, will 
imply an attraction mechanism towards General Relativity which starts to work during the 
radiation-dominated epoch. 



3. BBN in scalar-tensor theories 

We will now analyze the BBN processes, and the resulting primordial abundances, in 
the framework of scalar-tensor theories both with and without a self-interaction term. As 
stated above, we will consider throughout this section that the initial value of the scalar field 
is a free parameter, while its initial velocity is fixed to (pinit = 0. The initial values of the 
remaining variables are chosen so that they imply a flat Universe and lead, at the present 
temperature Tq = 2.725, to cosmological parameters given by: fi^.o = 0.27, ^a.o = 0.73 and 
Ho = 71 km s-iMpc"\ We assume an standard particle content of the Universe, with three 
families of light neutrinos, and the WMAP estimate of the present baryon to photon ratio 
(r] X 10^° = 6.14 ±0.25). 

The numerical computation consists of two main parts: first, we use a sixth-order Runge- 
Kutta integrator to fully determine the cosmological model and, in particular, the evolution 
of the expansion rate. Once determined the cosmological evolution in a given gravitational 
theory, we compute the BBN processes and the resulting primordial abundances of light 
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elements than ks to a complex networ k of 28 nuclear reactions, and usi ng the Beaudet and 
Yahil scher ae (IBeaudet &: Yahill 119771 ) . We used the reaction rates of ICaughlan fc Fowler 
(I19981) and ISmith &: all (119931 ) and the updated reaction rates of the NACRE collaboration 



( lAngulo &: allll999l ). The only aspect that differs from the standard BEN in General Rela- 



tivity is the expansion rate of the Universe. It impacts on the various nuclear abundances 
because they are determined by many reactions whose efficiency depends on the ratio of the 
reaction rates Fj and the expansion rate H. Indeed, a reaction with rate Fj is in equilibrium 
when Ti/H > 1, whereas the reaction is frozen when Ti/H < 1 . This is directly linked 
with the fact that the expansion dilutes the particles. Then, a modification in the expansion 
history may significantly change the nuclear abundances by modifying the dynamical struc- 
ture of the network of reactions, making some reactions more efficient, and limiting others 
(it should be noted that this remark is valid for the nuclear reaction rates, and also for the 
rates of the weak interaction processes that take place before BEN and that interconvert 
neutrons and protons). 

In order to characterize the deviation from General Relativity, it is convenient to in- 
troduce the speed-up factor, defined as the ratio between the expansion rate H{T) and the 
corresponding expansion rate in General Relativity Hgr{T) at the same temperature: 

H{T) 



(17) 



HcniT) 

When ^(T) > 1 (^(T) < 1 ), the Universe expands faster (slower) than in General Relativity 
at the temperature T. 



3.1. Theories without a self-interaction term 

We will first consider the simplest case of scalar-tensor models without a self-interaction 
term. In this case, the effective potential during the radiation-dominated epoch is (see 
equation f|T3l) ): 

K//(^)=0 (18) 

and, since we are assuming the initial condition f|T6l) . the scalar field will be frozen to its initial 
value (pinit (except for a slight temporary perturbation during the annihilation of electrons 
and positrons around T = 0.1 MeV, when Wb slightly deviates from 1/3). Consequently, 
= during all the EEN processes, and the relation flTTl) reduces to H = H^/A{ipi). 
Introducing the expression for in this equation, we can write: 



= V 3 P 



(19) 
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which imphes ^ ~ ^if i) = constant. 

Equation (fT9l) shows that the dynamics of the Universe is the same as in General Rel- 
ativity, but with an effective gravitational constant given by Gg// = G^A \ where = On 
is the usual Newtonian value {Gm = 6.673 x 10~^^m^ kg~^ s~^ ). Cyburt (jCvburtll2004l ) has 
determined the scaling of the various primordial abundances in terms of the physical con- 
stants. In particular, when the dynamics of the Universe is governed by General Relativity, 
the scaling with the effective gravitational constant is: 

0.35 

= 0.2484 ( ) (20) 



G 



N 



D/H 
^He/H 



7 



Li/H 



2.75 



8.65 



3.82 



eff 



Gn 
Gn 



0.95 



10" 



0.34 



10" 



G 



N 



'0.72 



10 



-10 



(21) 
(22) 
(23) 



Figure [T] shows these abundances as functions of Geff/G^, as well as their observational 
constraints. One clearly notes from this figure that there is no consistent value for Geff/GN 
that can simultaneously account for all the observed abundances: the observed ^Li primordial 
abundance requires G^ff/GN > 1-8 while ^He and D impose G^ff/GN ~ 1- 



1 1 1 1 1 1 1 1 1 


15 - 
10 - 


^ : 

.1.1.1.1. 


5 - 





Fig. 1. — Primordial abundances as functions of Geff/GN- 
assumed observed abundances. 



The shaded regions are the 



Consequently, in absence of a self-interaction term, scalar-tensor theories with an arbi- 
trary coupling function a{ip) (with (pinu = 0) cannot constitute a solution for the lithium 
problem. The same conclu sion can be found using the formulation of Kneller and Steigman 
feneller fc Steigman|[20oi ) 
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3.2. Theories with a self-interaction term 

We will now consider scalar-tensor models with a self-interaction term. As emphasized 
above, such a term acts as an effective potential even during the radiation-dominated epoch. 
Therefore, except for very early times (with a very high density), the scalar field does not 
remain frozen to its initial value and the speed-up factor is not necessarily a constant during 
the BBN processes. 

Of course, depending on the functional form of a{Lp) and V{lp), there exist many theories 
that exhibit a varying speed-up factor. Before analyzing a particular choice for a{(f) and 
V{(f), we may ask whether we can constrain the behavior of the speed-up factor in order to 
solve the ^Li problem without affecting the other abundances. 



3.2.1. Computation of the primordial abundances 

Since essentially all the neutrons are incorporated into ^He nuclei, the final production 
of ^He is directly related to the abundance of neutrons in epochs prior to the BBN processes. 
Indeed, Yp is roughly given by 

^ J^Mf^ (24) 

where (n/p) fr denotes the neutron to proton ratio when their weak reactions, n+Ug <-> p+e~ 
and n + ^ p + Ug, freeze out of equilibrium. Such a ratio is then fixed at a temperature 
Tfr of about 1 MeV, well before the synthesis of light elements (from ~ 100 KeV to ~ 10 
KeV). The {n/p)fr ratio, and hence the final Yp value, strongly depends on Tfj. through: 

{n/p)fr^eM-Q/Tfr) (25) 

where Q is the difference between the neutron and proton masses. A faster (slower) expansion 
rate of the Universe implies a higher (lower) freezing-out temperature, resulting in higher 
(lower) neutron and ^He abundances. 

The ''Li production is instead determined by the processes of BBN. In particular, when 
the baryon to photon ratio is greater than 3 x 10"^'^, the lithium is mainly produced through 
the reaction '^He(a,7)^Be followed by a decay of ''Be through electron captures (when ri < 
3 X 10^^*^, the lithium is instead produced by ^He(Q;, 7))^Li). Therefore, the final abundance 
of lithium strongly depends on both the helium production and the efficiency of the ^He 
burning to give ^Li. If Fj denotes the reaction rate, such an efficiency is given by Ti/H — 
{Ti/HRc)/i{T), so that a faster (slower) expansion rate of the Universe during BBN imphes 
a less (more) efficient production of ^Li. 
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Since General Relativity leads to a predicted Yp value in agreement with observations, 
one can expect that any other gravity theory must imply a speed-up factor close to unity at 
T ~ 1 MeV, or slightly smaller than unity to favor a less efficient production of ^Li without 
implying an over-production of ^He. Later, during the BBN processes, a speed-up factor that 
has significantly increased above unity could help to solve the lithium problem. Obviously, 
since any gravity theory must converge at late times towards General Relativity (in order to 
be compatible with Solar System experiments), the increase of the speed-up factor must stop 
at some time and, afterwards, it must approach unity. Therefore, models implying a constant 
or monotonic ^{T) function are not good candidates to solve the lithium problem. We will 
then restrict our analysis to self-interacting scalar-tensor theories satisfying the following 
two additional conditions: 1) they imply a speed-up factor with a local maximum (it is 
a non-monotonic function of T), 2) they have an attraction mechanism towards General 
Relativity. 

The analysis of the early behavior of £ is a complex problem, spe cially when (p starts 



to deviate from zero (see, e.g., (lAlimi fc Senalll997l : iNavarro fc all 120021 ) for theories without 



a self- interaction term). Nevertheless, the first condition (non-monotonic ^) is more easily 
satisfied in theories implying ^ < 1 at very early times (when ipmit ~ and the self interaction 
term is still negligible). Since C, ~ A{{p) at such early times, the condition < 1 implies 
A{ip)init < 1- Taking into account equation ([3]), the simple choice of a power law a{(p) oc y^^" 
for the coupling function, will imply ^inu < 1 provided that ipinu < 0. On the other hand, as 
quoted in section a similar choice V{(p) oc ip"^"^ for the self interaction term, will imply the 
existence during the radiation-dominated epoch of an attraction mechanism towards General 
Relativity (^ will not remain frozen to its initial value). Nevertheless, a self-interaction 
corresponding to m = 1, in other words, a simple mass term, is not a viable choice because 
the scalar field is not damped efficiently. In fact, by a simple dimensional analysis of equation 
(IT^ . one can define a characteristic time for the friction: Tf^ic = 'Irjj-i a characteristic 

' m, 

'^drag 

same order of magnitude (which is required by the fact that the self-interaction must play a 
role at temperatures relevant for BBN), one finds, when (y9 -C 1 that: 

Tdrag \ ^ j Tdrag 



time for the dragging in the effective potential: T^rag = \J jFJf\- When comparing the ratio 
for a model with m > 2 and a model with m = 1, if the constant coefficients are of the 



'^fric J m>2 V '^fric J ^=1 

That means that the friction is less efficient in the case m = 1, when the scalar field has 
converged in a neighbourhood of 0. This is supported by numerical computations that 
show that the scalar field quickly reach a neighbourhood of 0, but is not damped enough 
during the radiation dominated era with a self-interaction of the type ^(v?) oc pP' . In such a 
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case, the energy contribution of the scalar field becomes more and more important until the 
beginning of the matter dominated era. Such models might eventually lead to a very different 
estimation of the baryon to photon ratio through CMB anisotropies, then invalidating the 
hypothesises at the basis of the computations presented in this paper. 

In order to analyze the possible implications of scalar-tensor theories with a self-interaction 
term, we will adopt here the simple choice: 

a{^) = aip^; V{ip) = A^^ (26) 

with i^init = and ipinu < 0. 

The effective potential for these theories writes: 

W..A)^(l-3..)^,' + jj!^. (27) 

which means that, during the radiation dominated period, Veff{(p,X) = A'^ip'^/ (AnG^p^), 
leading to an attraction towards that is effective as soon as A^y9^ ~ 47rG'*p*(A), or in other 
words, as soon as the energy density of the scalar field is of the same order as the energy 
density of radiation. 

Using fl26|) . and varying the parameters {a,A,ipinit), we have numerically computed the 
BBN processes and the resulting primordial abundances. We have found that all models with 
— 1.5 < (fiinit < —0.9 lead to primordial abundances that agree with all the observational 
constraints mentioned in the Introduction. Consequently, there is not a lithium problem for 
such models. The parameters a and A are not strongly constrained. The only requirement 
is that A should not be too small in order for the attraction mechanism described above 
to occur during BBN, at temperatures between 10 MeV and 1 MeV. Indeed, from equation 
f ll2p . we know that the attraction mechanism approximately (when \(p\ ~ 1) occurs when 
A^ ~ G'*p^,(A); for 1 MeV < T < 10 MeV, which correspond to the period of BBN, since 
p*(A) oc (T/To)^ we are left with ~ 0.1 s^^ < A <~ 6 

To be more precise, we now present a special case by fixing the initial conditions to 
0mit = and ipinit = —1.3. Figure displays, in the (a. A) plane, the regions of theories 
leading to the observed primordial abundances. The constraints only come from the ^He and 
D abundances, which means that the ''Li abundance is in perfect agreement with the obser- 
vations for a much wider range of parameters. The two se parated regions corre s pond to the 



two different constraints imposed on the He abundance (iLuridiana fc all 12003 



Izotov fc al 



19991). 



The dynamics of the scalar field and of the speed-up factor are respectively shown 
in figures ([3]) and (jl]) for the particular choice (a. A) = (1,0.3 s^^). In a first phase, the 
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a 



0.1 0.1 0.J 0.4 0,5 0.6 
A 



Fig. 2. — Space of the parameters (a, A) for the theory defined by a{ip) = aip'^ and V{ip) = 
A'^if^ with fixed initial conditions: ipi = —1.3 and ipi = 0. The acceptable theories are 
represented by the shaded regions. 

self-interaction plays no role and the field is almost constant until T ~ 5 MeV, leading 
to an expansion slower than in General Relativity. Then, the y?'^ term starts to dominate 
the effective potential and attracts the scalar field towards zero. The expansion factor also 
oscillates and is driven to values greater than those obtained in General Relativity. The 
primordial abundances predicted in this case are Yp = 0.2449, D/H = 3.030 x 10~^ and 
^Li/H = 2.387 x 10^^'^, in agreement with observations. 

The clues for the success of this kind of scalar-tensor theories in solving the lithium 
problem have been already advanced at the beginning of this subsection. The evolution 
of the (n/p) ratio, as well as that of the ^He overproduction in the scalar-tensor model 
exemplified here, are shown in figure ([5]). As can be noticed on that figure, the effect of the 
speed-up factor on the n/p ratio is an integrated effect at temperatures between a few MeV 
and a few hundreds of keV. The fact that the speed-up factor is less than 1 at temperatures 
higher than 1 MeV implies that the n/p ratio departs from thermal equilibrium slighly later 
than in General Relativity, and this effect is partly compensated by a speed-up factor greater 
than one after 1 MeV, resulting in a n/p ratio that decreases faster in General Relativity 
than in the scalar-tensor model between 1 MeV and 0.1 MeV. These two effects almost 
compensate each other and result in a n/p ratio at the beginning of BBN comparable in the 
two models, and so, in a ^He abundance in the scalar-tensor model that is not very different 
from the one obtained in General Relativity. During the whole BBN period, the speed-up 
factor is instead significantly higher than unity. Then, the burning of ^He to produce heavier 
elements is less efficient than in General Relativity, and the resulting final ^Li abundance 
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Fig. 3. — Evolution of the scalar field ip during BBN for the scalar tensor theory defined by 
a{(p) — (/?^ and V{(p) — 0.3^</?^ with (pi — —1.3 and cpi — 0. 

( 

2-1 1 





Fig. 4. — Speed-up factor as a function of the temperature for the scalar tensor theory 
defined by a{(p) — </7^ and V{(p) — 0.3^^3^ with cpi — —1.3 and (pi = 0. The speed-up factor 
converges towards 1 as soon as the matter-radiation equality is reached. 
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is low enough to be consistent with the observed value. Indeed, the reaction ^He(a, 7)'''Be 
dominates the creation of ''Be (and then of ''Li by /3-decay of ^Be) during BBN. Therefore, 
neglecting the destruction of ^Be through the reaction ^Be(n,p) ''Li that is inefficient because 
of the low neutron abundance, a rough estimate of the final abundance of these elements can 
be obtained from dYser/dT = TYHe^yHeA/iTH) or, for small variations: 



where V is the reaction rate of ^He(Q;, 7)^Be, while iBe?, YHez and Yhca denote the abundances 
of ''Be, ^He and ^He, respectively. 

The lower panel of figure E] shows the evolution of the ''Be+^Li abundance in this scalar- 
tensor model as well as in General Relativity (GR). It can be seen from this figure that both 
predictions start to significantly deviate only at temperatures smaller than 0.06 MeV, and 
reach their final constant values at about T ~ 0.03 MeV. By applying the equation (1281) 
within such a small range of temperatures, we have: 



that also provides an estimate of the final Yu-jjYj^^ ratio. At temperatures below 0.06 MeV, 
the abundances of ^He and ^He are constant in both models, with Y^es ~ ^^"^ ^He^ ~ 
y^g4, while the speed-up factor oscillates around ^ ~ 1.16. Therefore, equation fl29|) implies 
^Lij ~ 0-8^Li7^; in rough agreement with our result obtained from the full integration of the 
nuclear reaction network: Yui = 0.74y^^. 

Of course, other theories with a different coupling function and/or a different potential 
could have the same effect on the ^Li abundance, provided that they lead to a speed-up factor 
with a time evolution similar to that found in the models analyzed in this section. As an 
example, the Appendix shows that, when the ipmit = condition is relaxed in scalar tensor 
theories without a self-interaction term (then, the convergence towards General Relativity 
is not ensured, and requires additional conditions), some of the resulting models can also 
solve the lithium problem. Therefore, the mechanism found here, based on a specific shape 
of the variation of the expansion rate during BBN, can be achieved in a great variety of 
scalar-tensor cosmological models. 



AFBe7= T7^He3lHe4AlnT 



(28) 




(29) 




Fig. 5. — Upper panel: Evolution of (n/p) as a function of the temperature in General 
Relativity (dashed line) and in the same model as in Figure H] (solid line). The dashed 
dotted line represents the (n/p) ratio at equilibrium. Medium panel: Evolution of the 
difference between "^He abundances in General Relativity and in the same model as in Figure 
m Lower panel: Evolution of the abundance of ^Be+^Li as a function of temperature in 
General Relativity (dashed line) and in the same model as in Figure H] (solid line). 
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3.2.2. 



Beyond BBN processes: effects on CMB and matter power spectrum 



It is important to note that, in the analysis of the previous subsection, at the end of 
BBN, a{ip) < 10~^ and A{Lp) ~ 1, so that Geff/GN ~ 1; moreover, a{ip) goes on decreasing 
because ip decreases, so that one can consider that gravitation is indistinghishable from 
General Relativity already at the end of BBN. Nevertheless, the scalar field still has a 
non negligible energy density, and consequently, the expansion rate is not the one of General 
Relativity with matter dust and radiation fluids, but receives an additional contribution from 
the scalar field energy density. This is why the speed-up factor is not immediately 1 but 
oscillates around a value of 1.16 until the end of the radiation dominated era and converges 
towards 1 when dust matter begins to dominate. As a consequence, we expect that such a 
model will have some impact on observables such as the matter power spectrum or the CMB. 
This is actually the case, but, these impacts are really of a different n ature that the ones that 
were previously investigated in the context of scalar-tensor theorie s (jChen &: Kamionkowski 
I999I : iBaccigalupi fc alll2"oOol : iTorreslbood : iNagata fc aJbood boosh : these works address the 
effects of modification of gravity (mainly variations of the gravitational coupling) on the 
CMB and matter power spectrum. In the models presented in this paper, the variation of 
the gravitational coupling plays a very important role during BBN (as examplified above), 
well before the matter dominated era, but at the epoch of matter-radiation equality, the 
gravitational coupling no longer significantly varies, and the scalar field behaves almost like 
a standard scalar field (a(v?) ~ 0). Let's know consider the impact of our model on the 
matter power spectrum and on the CMB. The matter power spectrum turn-over is defined 
by the scale ente ring the Hubble horizon at the matter radiation equality (that is H~^) 
(ICoble fc allll997l ). Since the expansion rate is slightly faster in the model presented here 
tha n in General Relativity at the time of equality, we expect a shift in the turn over given 
by JBacciealupi fc allboooh : 

eg 



Which yields, if we define the speed-up factor at equality by ^eq- 



6hi 
kr 



U - 1. 



(31) 



With E,eq ~ 1.075, one then has: Skr/k r ^ 7.5%. This sma ll shift is completely compatible 



with the current incertainties on data (ITegmark fc all 120061 ). 



The position of the first acoustic peak in the CMB power spectrum is also affected. 
The acoustic oscillations occur at an angular scale proportional to the size of the CMB 
sound horizon at decoupling and inversely proportional to the comoving distance between 
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the observer and the last scattering surface. So, if one defines the conformal time: 

dt 

Jo a{V 

the multipole moment associated with the first peak is given by: 

; T"o Tdec /oo\ 
Ip OC [66) 

To 

where tq is the conformal time of the observer (today), and r^ec the conformal time at 
decoupling. Then the shift in the position of the first peak: 6lp/lp = {Ip — lp^)/lp^ can 
be numerically inferred from the simulation. For the model p resented above, one finds: 



5lp/lp ~ 8%. This is a rather important shift (iPage fc alll2003l ). but it should be stressed 
that it depends strongly on the value of the expansion rate today, and that it has been 
assumed, to make a significant comparison that this expansion rate today was the same in 
the two models (General Relativistic cosmology and scalar-tensor cosmology). It is possible 
to play with the acceptable value of i^o (that is to give different values of Hq to the two 
models) to make this shift smaller. 

Finally, one expects that the model presented above will lead to other distortions in 
the CMB power spectrum at small scales since the expansion rate is different from the 
one of General Relativity during the whole radiation dominated era, when the small scales 
(those smaller than the one of the firsts acoustic peak) are entering the sound horizon, but 
investigating these problem demands a detailed analysis of the CMB physics in this new 
context that is beyond the scope of this paper, and will be dealt with in a forthcoming 
paper. 



4. Conclusion 

In the framework of the standard BEN, the WMAP measurement of the baryon density 
of the Universe Vt^h? = 0.0224 ± 0.0009 leads to a serious discrepancy between the predicted 
and observed ^Li abundance, even when systematic errors are cautiously taken into account. 
We addressed this problem by renewing the standard BEN scenario in a simple way: by 
considering that gravitation is described by a scalar-tensor theory. The expansion of the 
Universe during EEN is then modified, then modifying the conditions and the efficiency of 
nuclear reactions while assuming the standard nuclear physics. We showed that it is possible 
to obtain a small enough ^Li abundance and, at the same time, to preserve the ^He and 
D abundances thanks to a scalar-tensor theory of gravity with a non monotonic speed-up 
factor that has a generic behavior. The expansion must be similar to that found in General 
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Relativity at T ~ 1 — 2MeV, so that the freezing-out temperature of the weak processes 
n + ^ p + e~ and n + ^ p + Ue is also similar to that obtained in General Relativity, 
then leading to almost the same ^He production. However, during BBN, the expansion must 
be faster than in General Relativity, making the reaction ^He(a,7)^Be less efficient to burn 
'^He and to produce ''Li. 

Consequently, the ^Li abundance obtained for the baryon density inferred from WMAP 
could be a good imprint for the presence of a scalar field explicitly coupled to matter at the 
beginning of the Universe, challenging our understanding of the primordial Universe and the 
nature of gravity at early epochs. It could then be valuable to study the possible signatures 
of this scalar field in infiationary scenarios and to perform a detailled analysis of its impact 
on CMB physics and structure formation. 
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APPENDIX: Varying speed-up factor in theories vi^ithout a self-interaction 

term 

All the models analyzed in this paper starts from (pinu = 0. As a consequence, in any 
theory without a self-interaction term, the scalar field has a constant value during the whole 



radiation-dominated period. Nevertheless, previous works (lAlimi fc Senalll997l : iNavarro &: al 



20021 ) have shown that, when the condition (finit = is relaxed, many of these theories also 
imply a non-monotonic evolution of the speed-up factor. The convergence of these theories 
towards General Relativity has instead important differences with respect to that found for 
the models analyzed in this paper. In non-monotonic scalar-tensor theories without a self- 
interaction term, the convergence towards General Relativity is not ensured independently 
of the initial conditions. 

One of the main conclusions of this paper is that the lithium problem can be solved in 
many scalar-tensor theories, defined by different coupling functions and/or different poten- 
tials, but implying an evolution of the speed-up factor with the same general properties as 
those found in Section 13. 2[ In this Appendix, we explore the impact on BBN of a theory 
defined, in the Dicke- Jordan frame, by 

3 + 2uj{<j)) = + b (34) 
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with (l)init > 1 (this is a necessary condition for the speed-up factor to be less than unity at 
the beginning). 



In the hmit close to general relativity 



1), the above theory leads (ISerna fc alll2002 



Barrow fc Parsons! 119971 ) to an Einstein frame formulation with a{(^) oc Iv^l^^^ and ^init < 0. 
The effective potential defined in (fT3ll then becomes 



K//(<^) = sgn((^) 



2(1 - 3wh 



I I ^ 



(35) 



where sgn(x) = —1, 0, 1 if x is negative, zero or positive, respectively. This effective potential 
does not have a minimum, but a stationary point at (f = 0. So, if the scalar field does not 
have a high enough initial velocity, it will roll down its effective potential during the matter- 
dominated era, pushing the theory far from General Relativity. Then, in these theories, one 
must give an initial velocity to the scalar field so that i t can reach the oppos i te sid e of the 
stationary point before the matter dominated era. In iDamour fc Nordvedtl (119931 ). it was 
shown that, during the radiation dominated era, the total displacement of the scalar fields 
is given by: 



^ 2 



(36) 



1 - ^yv3^ 

So that we can deduce the initial velocity necessary to reach 0, starting from ipi. It is given 
by: 



1 + expi-2)ip,/V3 



(37) 



that is always finite and bounded by — and y^- 



Since it is not possible to perform the analytical transformation (jl]) to find its exact form 
in the Einstein frame, the models defined by (IMI) have been integrated in the Dicke- Jordan 
frame. Moreover, the integration was performed backward in time, then replacing the initial 
condition on the scalar field by the present value of the coupling function co'(0). The present 
value of the derivative of the scalar field was fixed to zero. 

By varying the three remaining parameters (a, 6, a;(0o)); we identified the cosmo logical 
models that lead to the observed abundances of light elements. The results were not sensitive 
to a;((/)o) as long as a;((/)o) is sufficiently large to pass the Solar System tests. Then, restrict- 
ing the analysis to the parameters a and b, the figure ([6]) shows (as shaded regions) the 
(a, b) couples leading to primordial abundances of ^He, D and ^Li compatible with observa- 
tions. Again, the two distinct adm issible regions a r e rela t ed to the two differ ent observations 
considered for the '^He abundance iLuridiana fc all (120031 ): llzotov fc all (Il999l ). 
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Fig. 6. — Acceptable theories are represented by the shaded regions in the (a, h) plane, for 
3 + 2u;(0) = + 6 in the Dicke- Jordan frame. 

Speed-up factor 




T (MeV) 

Fig. 7. — Speed-up factor as a function of the temperature (in MeV) during BBN. 

As an illustration of these theories, the figure ([7]) shows the evolution of speed-up factor 
corresponding to the particular choice a = 4, 6 = 2.68 and uj{(f)Q) = 10^^. It exhibits the 
behavior we described previously in order to solve the lithium problem and to preserve 
the ^He abundance. The predicted primordial abundances in this example were in fact: 
Yp = 0.2405, D/H = 2.871 x lO'^ and ^Li/H = 2.601 x 10"^°. 
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